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Abstract: Scalar multidimensionally consistent quadrilateral lattice equations 
are studied. We explore a confluence between the superposition principle for 
solutions related by the Backlund transformation, and the method of solving 
a Riccati map by exploiting two known particular solutions. This leads to an 
expression for the iV-soliton-type solutions of a generic equation within this 
class. As a particular instance we give an explicit 7V-soliton solution for the 
primary model, which is Adler's lattice equation (or Q4). 



1. Introduction 

Hirota's ./V-soliton solution of the Korteweg-de Vries (KdV) equation [1,2] was 
characterised in terms of a Backlund transformation originally by Wahlquist and 
Estabrook in [3] . Nowadays the close relationship between soliton solutions and 
Backlund transformations is a part of the landscape of soliton theory. 

The discrete systems which form the natural analogues of integrable partial 
differential equations of KdV-type are the scalar multidimensionally consistent 
quadrilateral lattice equations. Amongst the features which are more transpar- 
ent on the level of the discrete equations is the Backlund transformation, in fact 
a natural auto-Backlund transformation is inherent in the discrete equation it- 
self because of the multidimensional consistency. This is helpful in defining the 
very notion of soliton solution in the discrete setting, where the usual physical 
properties attributed to these solutions are less well studied. It also makes the 
solutions arising from iterative application of the Backlund transformation very 
well suited to a treatment on the discrete level. Particularly (we have found) 
when combined with an elementary technique for solving the Riccati equation, 
which we now take a moment to explain. 
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In the scalar case the discrete Riccati equation (Riccati map) is an equation 
of the form 

lZ(u, u) := ciq + a\u + a2?i + a^uu = 0, (1-1) 

where u = u(n) and u = u(n + 1) are values of the dependent variable u as a 
function of the independent variable n € Z, and in general the equation is non- 
autonomous so that the coefficients a . . . a 3 G C are functions of n. One of the 
many interesting properties of this equation is as follows: given two particular 
solutions v = v(n) and w = w(n), the substitution 



V — TjW 



(1.2) 



reduces (1.1) to a homogeneous linear equation for the new variable r\ = 77(71). 
Specifically 

lZ(w, v) _ . ,„ „. 

V = -^V => *(«,«) =0 (1.3) 

as a consequence of the vanishing of the first and last terms in the expansion 

(1 — 77) (1 — rj)lZ(u, u) = 1Z(v, v) — r]lZ(w, v) — rj1Z(v, w) + nfjlZ{w, w). (1.4) 

This technique has been applied to the Riccati map associated with the lattice 
Backlund transformation previously [4,5]. In the present communication we show 
that a similar technique applied to the superposition principle makes it possible 
to construct a closed-form expression for the solution found by the application 
of N Backlund transformations. In other words, it leads to an expression for the 
iV-soliton-type solutions. 

From classification results due to Adler, Bobenko and Suris (ABS) [6,7] there 
has emerged a core list of canonical forms for the scalar multidimcnsionally 
consistent quadrilateral lattice equations. With one notable exception, explicit 
iV-soliton solutions have been found for all of the equations on the ABS list [8, 
5,9], and their relation to the Backlund transformation has also been examined 
[5]. The exception is the lattice equation discovered by Adler [10] (cf. also [11, 
12]), which acquired the name QA in [6]. Q4 holds the position as the primary 
model on the ABS list because all of the others can be found from this single 
equation by a process of degeneration. The principal application of the technique 
developed here is therefore to Q4, for which an explicit A^-soliton solution has 
not (to our knowledge) hitherto been given. 

We start in Section 2 by specifying the class of equations under consideration. 
In Section 3 the Backlund transformation and superposition principle of a generic 
equation within this class are written in terms of the associated variable 77. The 
iterative application of the Backlund transformation in this new variable then 
leads to the main iV-soliton formula in Section 4. In Section 5 we provide a 
simple illustrative example, using equation Al from [6]. In Section 6 we use the 
A-soliton formula to exhibit an explicit iV-soliton solution of Q4, this builds on 
the seed and one-soliton solution constructed for this equation previously in [4] . 
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2. The class of equations 

We will study an equation of a generic form which encompasses all of those listed 
by ABS in [6]. Specifically, we consider a polynomial of degree one in four scalar 
variables, 

Qp, q (u, u, u,u) := ao(p, q) + a\(p, q)u + a 2 (p, q)u + ... + ai 5 (p, q)uuuu, (2.1) 

in which the coefficients are functions of lattice parameters p and q. We assume 
the Kleinian symmetry and covariance of the polynomial, 

Qp,q(u, U, U, ll) = Qp,q(u, U, U, U) = Qp.q(u, U, U, U), (2.2) 

Qp tq (u,u,u,u) = -Q qt p(u,u,u,u), (2.3) 
as well as compatibility of the system 

Qp,q(u, U, U, u) = 0, Qp,q(u, U, U, u) = 0, 

Q q ,r{u, U, U, u) = 0, Qq, r (u, U, U, u) = 0, (2.4) 
Qr,p{u, U, U, u) — 0, Qr,p(u, U, U, V,) — 0, 

for any p, q and r. By the compatibility of (2.4) we mean that given u, u, u, u, 
u, u and u satisfying the equations on the left, the remaining equations on the 
right are all satisfied by the same value of u. It is from the polynomial Q that 
we then define an autonomous quadrilateral lattice equation, 

Qp. q (u,u,u,u) = 0, (2.5) 

where u = u(n, to), il = u(n+ 1, to), u = u(n, m + 1) and u = u(n+ 1, m+ 1) are 
now values of the dependent variable u as a function of independent variables 
n, m G Z. The covariance (2.3) and the compatibility of (2.4) make the lattice 
equation multidimensionally consistent [13, 14]. 1 



3. Backlund transformation and superposition principle 

We are interested in the solutions of (2.5) constructed by the application of 
N consecutive Backlund transformations to an initial seed solution. To apply 
the Riccati solution technique we will presuppose that the seed solution admits 
covariant extension. 

Definition 1. A covariantly extended solution of (2.5) is one which also satisfies 
the (compatible) equations 

Qp,ii(u, u, Tjit, Tiu) = 0, Q q ,ii(u,u, TiU, Tjti) = 0, (3.1a) 
Qh,ij (u, TiU, Tju, TfTju) = 0, (3.1b) 

for all i,j e{l...N}. 

1 It should be expected that all polynomials satisfying (2.1) and (2.2) admit paramctcrisation 
leading to the multidimensional consistency [7,15]. On the other hand, for polynomials with 
no non-linear terms the assumption (2.2) seems less natural [16]. 
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In this definition we have introduced shift operators Ti . . . Tjv associated with 
new lattice parameters l\ . . . In, so the requirement is that u inhabits an (N + 2)- 
dimcnsional lattice satisfying the equations (2.5) and (3.1) throughout. We re- 
mark that situations in which there is no natural covariant extension are con- 
ceivable. However, covariant-extendibility is not a restriction because generically, 
given a solution of (2.5), the multidimensional consistency implies that solutions 
of the equations (3.1) exist. 2 

We give a definition of the (extended) Backlund transformation of (2.5). 

Definition 2. When two functions u and v satisfy the system of equations 

Q Pt r(u,u,v,v) = 0, Q q , r (u, u,v,v) = 0, (3.2a) 
Q u , r (u,T l u,v 1 T t v) = 0, ie{l...N}, (3.2b) 

we say they are related by the Backlund transformation of (2.5) with Backlund 
parameter r, and write u ~ v. 

One usually considers only the two-dimensional system (3.2a) in the definition 
of the Backlund transformation, but here we have included auxiliary equations 
so that the Backlund transformation preserves the covariant extension. The fol- 
lowing statement is locally equivalent to the assumption of multidimensional 
consistency described in Section 2. 

Lemma 1. If u is a covariantly extended solution of (2.5) then considered as a 
system of equations for v the Backlund equations u ~ v are compatible. Moreover 
the function v which emerges is again a covariantly extended solution of (2.5). 

Now, the special choice of Backlund parameter r € {li . . . In} causes an interplay 
between the Backlund transformation and the covariant extension, and it is this 
which facilitates the Riccati solution technique. 

Lemma 2. Let u be a covariantly extended solution of (2.5). Then for each 
i € {1 . . . N} the substitution 

... .PT 1 -^,]* (33) 



1 - Vi 

reduces the Backlund equations u ~ m to the following compatible homogeneous 
linear system for the new variable r\i, 

Qp,i i {u,u,T i u,J lU ) Q^KMjj u,TiU) 

Qi , iAu, TjU, "Hit, T^TV -1 ^) r , 

T ^ = ~ n T T-1 T T 3 Z{1 ■•■!*}■ 3.4b 

Qi 3 ,u( u > T 3 u ' J t u,TjTiu) 



2 All solutions of (2.5) arising as ^-cycles of the Backlund transformation - which are a 
discrete analogue of finite-gap solutions [17-19], do have a natural covariant extension (cf. the 
example in Section 6). On the other hand, for solutions which result from a generic initial 
value problem on the lattice, a natural covariant extension is not manifest. 
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Proof. The covariant extension provides two particular solutions of the system 
u Ui, namely Uj = TiU immediately from (3.1), and Ui — T~ 1 u, which can be 
seen by applying T^ 1 to (3.1) and using the symmetry (2.2). Now observe that, 
considered as equations for v, the Backhand transformation (3.2) is a system of 
compatible Riccati maps because of the assumption (2.1). The substitution (3.3) 

can then be seen to bring u ~ Ui to (3.4) by considering an expansion of the 
form (1.4) on each equation. The compatibility of (3.4) is inherited. □ 

We remark that the behaviour of rji and Ui away from the origin in the lattice 
direction associated with shift Tj is likely to be quite special. Actually in practice 
we will only use the covariant extension of u; t and rji into the remaining lattice 
directions associated with shifts Tj with j^i. 

Let us turn now to consider the superposition principle for solutions of (2.5) 
related by the Backlund transformation. Again we state the following without 
proof because it is a standard calculation based on the multidimensional consis- 
tency. 

Lemma 3. Let u be a covariantly extended solution of (2.5) and suppose u\ and 
U2 are such that u ~ u\ and u ~ u 2 for some n and r 2 . Then the function tti 2 
determined algebraically by the equation Q ri ,r 2 ( M : u ii u 2, "12) = satisfies the 
Backlund relations u\ ~ u 12 and u 2 ~ Ui 2 . 

This brings us to the crux of the method, it lies in considering Lemma 3 in the 
case that n, r 2 € {h ■ ■ ■ In}- The resulting interplay with the covariant extension 
(in which li . . . In play the role of lattice parameters, cf. Definition 1) leads to 
a relation, derived from the superposition formula, connecting now 77-variables 
associated with each Backlund transformation. 

Lemma 4. Let u be a covariantly extended solution of (2.5) and suppose the 
functions rji ... tjn are such that 

Ui = [T * 1 ~ ViJi]u u l rt Ui , ie{l...N}. (3.5) 

1 - Vi 

Then for each i,j£ {1 . . . N} the following identity holds 

(Tj^XTVfc) = (T-^KT^)- (3-6) 
Furthermore the new function rfij defined as 

Vij = ) ~ ^ TT'Vj (3-7) 



is such that 



[TJ 1 - rujlj\ui 
1 - mi 



(3.8) 



Proof. The proof is based on Lemma 3 plus a calculation which exploits the 
assumptions made about the equation (2.5) described in Section 2. First we 
observe that 



Q htlj («, Tf u, T>, Tf T» = 0, n,u e {+1, -1}, 



(3.9) 
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for all i,j € {1 . . . TV}. In the case \x = v = 1 this is immediate from the covariant 
extension (3.1b), the other cases follow from the symmetry (2.2). By combining 
(3.9) with the affine-linearity (2.1) we have for each /j,, v e {+1, —1} and each 
pair i, j e {1 . . . TV} 

Q iaj {u,ltu,l v jU ,z) = 5 y (u,T?u,T>)(z-T?T>) (3.10) 

as an identity in z. Here we have introduced a polynomial S of degree one in 
three variables defined in terms of Q as 

Sij(w,x,y) := Qi it ij{w,x,y, 1) - Q/^ (to, a;, y, 0) = d z Qi u i } {w,x,y,z). (3.11) 

This new polynomial inherits the symmetry 

Sij(w,x,y) = -Sji(w,y,x) (3.12) 

from the covariance of Q (2.3). The simplification which occurs by using (3.10) 
to replace occurrences of Q with S will be the principal mechanism exploited in 
the calculation. 

Now, according to Lemma 2 the assumption (3.5) means that each r\i satisfies 
(3.4). Using (3.10) in (3.4b) as well as in the equation found by applying T" 1 to 
(3.4b), we are able to deduce the following simplified equations for shifts on rji 
into the extended lattice directions, 

5 (,,^».T>) ^ *(«ffV«) (3 . 13) 
& i j{u,T i u,\jU) &ij(u,J i u, \ ■ u) 

which hold for all pairs i,j e {1...JV}. The relations (3.13) combined with 
(3.12) may be used to directly verify the property (3.6). 

We now proceed with the principal calculation; substituting for Ui and Uj 
leads to the following expansion of the superposition formula 

(1 - T)i)(l - Vj)Ql i ,l j ( U ' U i> U i> U ij) 

= Qi^i^uJ^uJ^u.u^) - r] i Qi ii i j (u,T i u,Jj 1 u,u ij ) 

- VjQh,i 3 ( u , ~T V Uij) + ViVjQhJi (u, TiU, Tju, mj), 

= S ij (u,Tr 1 u,Tj 1 u)(u i: j - TJ^TT 1 ^) - r] i S ij (u 7 T i u,Tj 1 u)(u ij ~ TjTT 1 ^) 

- r) j S ij (u,Tt 1 u,T j u)(u ij - T^Tju) + rji^jSij (u, Tjti, Tju)(uij - TiTju), 

(3.14) 

The first equality here results from substituting for and Uj using the equation 
on the left of (3.5) and exploiting the assumed affine-linearity of Q (2.1). To get 
the second equality from the first we have used (3.10). The third equality follows 
from the second by using (3.13) whilst bearing in mind (3.12). 
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The expansion (3.14) thus enables us to write an expression for the (covari- 
antly extended) solution of (2.5) determined from u, Ui and Uj by the superpo- 
sition formula Qi u i 5 (u, Uj, Uj,Uij) = 0, 

= [T^TJ 1 ( Tt^T.T- 1 - (J- 1 ^)!- 1 !, + (Tj ?fc ) (T ~ 1 t?.,- ) T t T 3 ] u 
1 - (T/V) - (T-\) + (TjVi^Vj) 

(3.15) 

It is easily verified that the expression for rjij in (3.7) leads to equality between 
given on the left of (3.8) and found by superposition in (3.15). The 
statement (3.8) therefore follows from Lemma 3. □ 

Thus it is demonstrated that, beyond what might be expected, the function rjij 
is determined completely in terms of r/j and rjj by an expression (3.7) which 
does not depend on the particular form of the polynomial Q. The additional 
statement of Lemma 4, (3.6), can be viewed as encoding the permutability of 
the Backlund transformations (specifically Ujj = Uji in (3.8)) as a property of 
the ^-variables; this property will also play a key role later. 



4. Backlund iteration 

Here we apply N consecutive Backlund transformations to a covariantly ex- 
tended seed solution of (2.5). Choosing the Backlund parameters to coincide 
with the lattice parameters l\ . . .1^ of the covariant extension (cf. Definition 1) 
renders the Backlund iteration scheme tractable and yields a closed-form ex- 
pression for the resulting iV-soliton-type solution. In the main theorem which 
follows it is both convenient and (as will be seen) natural to introduce a solution 
v and functions ...</> n which are shifted versions of the seed solution u and 
the variables r}\ . . . tjn of Section 3. 

Theorem. Given a covariantly extended solution of (2.5) denoted by u, let 

v=[I\l 1 T-i]u (4.1) 

and for each i e {1 . . . N} suppose the function <pi satisfies the system 

7 Qp,ii( T i v , T iV, T i v ,v) , 1 S g ^O,T t {;,Tft;,£) 

W = -7S Tt r -3 n^ft W = _ 7S rr =F^ TT^Vi, (4.2a) 

Q P ,h O, TiU, v, Tfv) Q q . h (Tiv, J l v, v, Tfv) 

Qij.M (Ti«, TiTjV, Tfv, Jjv) 
Qi^hO'i'", TiTjV, v, T?Tjv) 
Then the linear difference operators Bi . . . Bjv defined as 

B i: =l-&T?, ie{l...N}, (4.3) 

commute amongst each other, and the function defined in terms of these 
operators by the expression 



T;,o, „ ^ „, T2T ., fc. je{l...N}. (4.2b) 



"[nLB.]!' 



is a solution of (2.5) related to u by the composition of N Backlund transforma- 
tions with Backlund parameters l\. . .In- 
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Proof. Let us begin by supposing rji satisfies (3.4) for each i £ {1 . . . N}. Then 

according to Lemma 2 the function m defined by (3.3) satisfies u ~ Uj. In other 
words the following is true when we take i = 0: 

[T" 1 - Vi...ijT j ]u 1 ... i i . 

= — — — => ul.a ~ je{l...JV}. (4.5) 

t 

Now, using Lemma 4, if we define the functions recursively by the expres- 

sion 

,' Jf T, ',;,... ( ,; nj! i,je{l...JV}, (4.6) 

then (4.5) is true for any i = k £ {1 . . . TV} as a consequence of it being true 
when i = k — 1 . 

The recurrence (4.6) is not difficult to solve. Let us exhibit the first few 
iterations, writing them in the following way (which requires use of (3.6)) makes 
the emerging general formula quite evident: 

T j Pi -mTiJl 



„ _ T^MV^HV-mTi]! T -i T -i 

T 7 i T 2 -( T i ? 72)T 2 J[T 1 -jyiTiJl 



„ _ T,[T 3 ' 1 -(Tr 1 T 2 - 1 r ? 3)T3][T 2 ^-(Tr 1 r 72 )T 2 ][Tr 1 -r; 1 T 1 ]l x x x 
" 123 ^T 7 1 [T 3 - 1 -(Tr 1 T 2 - 1 ?73 )T 3 ][T 2 - 1 -(T r 1 ??2 )T 2 ][T r 1 -, 1 T 1 ]l l1 ' 2 ' 3 ^ 

To write the term in this sequence we introduce linear difference operators 
A -=il[ T 7 1 -([ni= 1 iT fe - 1 ]^)T,], i£{l...N}, (4.7) 

3 

in terms of which 

^■^[nuv]^. (4.8) 

This can be verified as the solution of (4.6) directly by substitution. 

We now reconstruct the solutions using the equation on the left of (4.5), in 
fact it is sufficient to assume i and j are consecutive, so we need only consider 
the recurrence 

M ,.,= [Trl ^ 1 "^ K -^- 1) ; i£{l...N}. (4.9) 
The solution is simply 

AiU 



Ml 



, i = _i_ ie {i...JV} (4.10) 



which, bearing in mind (4.8), can be easily verified by substitution. Thus we 
have found an expression for solutions u L A which satisfy the desired Backhand 

relations ~ tti...i, « € {1 . . . TV}. 
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Now, we began the proof by assuming that, for each ie {1 . . . N}, r\i satisfied 
(3.4). Actually if we make the association 

then direct substitution shows that (3.4) is just a consequence of (4.2) and (4.1) 
which occurred in the hypotheses of the theorem. Using the association (4.11) 
it is straightforward to see that 

A^tnLB.nntiTr 1 ]- (4-12) 

The commutativity of the operators Bi . . . B^v, which enables us to remove the 
ordering of their product in (4.12), follows from the identity 

&(T?&)=^(Ti&). (4-13) 

which is just (3.6) expressed in the shifted variables through (4.11). From (4.12) 
and (4.1) it is clear that in (4.4) is simply u\...n in (4.10) expressed in the 
shifted variables, which completes the proof. □ 

We conclude this section with a few remarks. The functions <j>i...<f>n arise 
in the commuting operators defined in (4.3) making it natural to adopt these 
variables in preference to r}\ . . . t]n. The function v in (4.1) is one of 2 N particular 
solutions of the iterative Backlund scheme which are provided by the covariant 
extension. The others are manifest in (4.4) when the functions (f>\ . . . (pN are 
taken to be the trivial solutions of (4.2), namely constant at or oo, implying 

U (N) g {[T -1 . . . j-l ]u = Vi [TlT -l . . . T -l] Uj . . . , [ Tl • • • J N ]u}. 

Note that the function in (4.4) lies on an (N + 2)-dimensional lattice, 
the desired iV-soliton-type solution as a function of n and m alone is realised 
by discarding the covariant extension, i.e., by evaluation at the origin in the 
additional lattice directions associated with the shifts Ti . . . Tjv- 

For the examples, as will be seen in the following two sections, it is useful to 
introduce the function 

/:=[Il£iBi]l, (4.14) 

which appears in the denominator of the iV-soliton formula (4.4). The operator 
present here has the expansion 

N N N N 

nBi = i-E*j?+E E *aton?Tl 

2—1 2—1 2—1 

N N N 

-EE E 0i(T^)(T^ fe )T^T 2 fe + ..., (4.15) 

i— 1 k—j+1 

so the function / is reminiscent of the one introduced by Hirota [1,2]. Addition- 
ally, we observe that the function / can be expressed in a factorised form 



/ = (1 - 0i..jv)(1 - 0i..jv-i) ■■■(!- 0i 2 )(l ~ <h) (4-16) 



10 



J Atkinson and F W Nijhoff 



through functions <j>i...i, which, for each i E {1 . . . TV}, are related to r\\...i present 
in the proof (4.8) by the relation 4>i...i = [IljLi+i ^J 1 ] Vi—i- O r m terms of the 
shifted variables alone 

, [T?n,-',Bj]i [n;-iBj]i 



5. An explicit iV-soliton solution of Al 

In this section we provide an illustrative example of the main theorem, in par- 
ticular the iV-soliton formula (4.4). The solution we consider was (up to a gauge 
transformation) discovered previously by other methods in [5]. 
Consider the expression 

_ „ ~ _ (u + u)(u + u) (u + u)(u + u) S 2 a i (p 2 - g 2 ) 

where a e C \ {0} and S e C are arbitrary constants. This is a polynomial of 
the form (2.1); direct calculation will verify that it satisfies all of the properties 
listed in Section 2. The equation defined by (5.1) is equivalent to Al in [6], 
the only difference being a point transformation of the lattice parameters p — ► 
a 2 /(p 2 — a 2 ), q — ► a 2 /(q 2 — a 2 ), which is useful for discussing solutions [5]. 
The equation defined by (5.1) has an elementary solution, u, of the form 

u = Aip + B/tp, (5.2) 

where 

and A, B G C are constants subject to the constraint 

AB = (S/4) 2 . (5.4) 

We will take u in (5.2) as a seed solution to start the Backlund chain. This 
solution admits a natural covariant extension; observe that 

7 a+ P I I a + ( l I IK K\ 

ip= V, ip= ^ (5-5) 

a — p a — q 

and consider complementing these with the compatible equations 

T^ = ^V, i€{l...N}, (5.6) 

(X Li 

in which we have introduced new parameters h . . . In € C \ {—a, a}. The equa- 
tions (5.6) extend ip to the (N + 2)-dimensional lattice (note that the function 
appearing in (5.3) should now be interpreted as the evaluation of tp at the origin 
of the extra directions). Of course the solution u defined in terms of ip by (5.2) 
also inhabits the (N + 2)-dimensional lattice, and it is straightforward to see 
it satisfies the additional system of equations (3.1). Thus, for the quadrilateral 
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lattice equation defined by the polynomial (5.1), we have given a solution (5.2) 
and its natural covariant extension through (5.6). 

We are interested in the result of N applications of the Backlund transfor- 
mation to the seed solution u. To apply the iV-soliton formula (4.4) we should 
construct from u the related solution v and the functions <pi . . . </>jv- The solution 
v is given in terms of u by (4.1). With u as in (5.2) we use (5.6) to find 

v=[Uli^] u = A'i, + B'/i; (5.7) 
where the new constants A' and B 1 are related to A and B by the equations 

N i N i ; 

i=l i=l 

but satisfy the same constraint (5.4). So u and v are the same up to a change in 
the value of some constants present in the solution. The other ingredient present 
in the TV-soliton formula (4.4) is the set of functions ...</> at, these are defined 
in terms of v by the system (4.2) for each i e {1 . . .N}. Direct substitution of 
(5.1) and (5.7) into (4.2) followed by the use of A'B' = (S/i) 2 , (5.5) and (5.6) 
yields the following system for each fa, 



P-k , i Q-h 



p+h q+h 



(5.9a) 



T i4>i = f—T<t>u je{l...N}. (5.9b) 

3 * i 

Integrating the equations (5.9a) we can write 

( f> i = < f> i (n,m) = <p i , (^f) (^) , ie{l...N}, (5.10) 

where <p\fi . . . 4>n,o are independent of n and m. 

So we have constructed the solution v and functions 4>\ ■ ■ ■ 4>n present in the 
iV-soliton formula (4.4). We proceed by using (5.6) and (5.9b) to write (4.4) 
in terms of un-shifted functions tp an d <f>i ■ ■ -<t>N\ first consider the function / 
defined in (4.14) which appears in the denominator of (4.4), we find 

N N N 

f = /(0 1: . . . , 4> N ) = 1 - & + E E 

N N N N N 

"E E E ^AkXlxlX^ + .^ + i-lfW^ J] Xfj), (5.11) 
where 

X ij := l f^i = -X ji (5.12) 

are constants such that Tj(f>i = Xij4>i. A similar consideration of the numerator 
in (4.4) yields, after collecting terms in tj) and 1/tp, 

=A'^ + (B'/^, (5.13) 
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where we have introduced the functions / which are defined in terms of / in 
(5.11) by the equations 

CI ~\~ li 

(5.14) 

So, through the functions /, /+ and / given in (5.11) and (5.14), given in 
(5.13) is a rational expression in the un-shiftcd functions ip and <pi . . .(f) m- Taking 
these functions as in (5.3) and (5.10) respectively, (5.13) thus gives the solution 
as a function of n and m alone. This is the 7V-soliton solution which results from 
N applications of the Backlund transformation to the seed solution u given in 
(5.2). It is worth observing that for (5.13) to be a true AT-soliton solution, i.e., 
one containing N constants of integration, it is necessary that Zj ^ for all 
ie{l... N} and that Zj ^ ±Zj for all distinct i,j e {1 . . . N} (these conditions 
are in addition to the one already mentioned, U {—a, a} for all 4 G {1 . . . N}, 
which was required for the covariant extension). 

Arranging the solution as we have done in (5.13) is suggestive of further struc- 
ture behind this expression. This further structure has a natural interpretation 
in the contexts of the Cauchy-matrix approach described in [5] (where this so- 
lution was first given) and the Casorati approach [8,9]. To finish this section we 
describe how this structure should be interpreted in relation to the approach in- 
troduced in Sections 3 and 4. The main idea is to extend the functions <p\ . . . 4>n 
into one further special lattice direction. This direction is associated with lattice- 
parameter a and we denote shifts in this direction by T , 

T a <Ai = —j-f <Pi> ie{l...N}, (5.15) 

which inspecting (5.14) enables us to write = Tj 2 /. (Note that to be consis- 
tent ip should take values in {0, oo} away from the origin in this special direction, 
but this won't play a role in the present construction.) The form of the solution 
(5.13) then motivates the introduction of a new linear difference operator 

X:=AVT- 2 + (B'/^T2 (5.16) 

for which the following facts can be verified 

XI = v, [Bj,X]=0, ie{l...N}. (5.17) 

The first is clear from the definition (5.16) whilst bearing in mind (5.7), the 
second follows from (5.16) and the definition of the operators Bi . . . B N in (4.4) 
because 

^(Tq 2 ^) = 4>i(Tfip), (l/VO(Tfo) = 0i(l/Tfy), ie{l...N}, (5.18) 

which are immediate from (5.15) and (5.6). Now, using the properties (5.17) in 
(4.4) and the definition (4.14) we find 

v.™ = ^. (5.19) 

Thus, for this example, have uncovered the presence of an operator X defined in 
(5.16) which, due to the properties (5.17), enables the iV-soliton formula (4.4) 
to be written in the simplified form (5.19). 



A constructive approach to the soliton solutions of lattice equations 



13 



6. An explicit 7V-soliton solution of QA 

We now use the 7V-soliton formula (4.4) to give an explicit TV-soliton solution 
for the equation listed as QA in [6]. The construction of a seed solution for QA 
is itself an interesting problem which was solved previously in [4], where the 
ensuing one-soliton solution was also given. The full iV-soliton solution which 
results from this seed solution is given here for the first time. 



6.1. The equation QA. We consider QA in the Jacobi form [20], which is related 
by a change of variables to the Weierstrass form given originally by Adlcr [10] 
(and expressed more concisely in [11]). It may be defined in terms of the poly- 
nomial 

<2p, q (u, it, it, u) := p(uu+uu)—q(uu+uu) — ^— — ^-y (uu + uu — pq(l + uuuu] 

(6.1) 

In Section 2 we did not specify the set from which the lattice parameters were 
taken, here p = (p, P) and q = (q, Q) are points on an elliptic curve, p, q e Y, 

r = r(fc) := {(x,X)\X 2 = l + x 4 - {k + l/k)x 2 } , (6.2) 

where k € C\ { — 1,0,1}, the Jacobi elliptic modulus, is a fixed constant. The 
properties listed in Section 2 may all be verified directly for the polynomial (6.1). 

In our consideration of solutions for QA a central role will be played by the 
natural product that turns Y into an abelian group. We recall some facts and 
establish some notation regarding this here. To start with let us write the rational 
representation of the group product: 

f pQ + gP Pp(g 4 -1)-Qg(p 4 -1) \ 
Pq il (1 - p 2 q 2 )(qP - pQ) I ' 



and note the group identity is the point e = (0, 1) and the inverse of a point p 
is p _1 = (— p, P). This group structure on Y is naturally parameterised through 
the Jacobi elliptic functions by introducing the mapping 

f: z (^fksn(z; k), cn(z; fc)dn(z; k)^j . (6.4) 

By convention (see for example Chapter 5 of [21]) the primitive periods of the 
function z t— ► sn(z; k) arc denoted AK and 2iK' (here i denotes the imaginary 
unit). The mapping f defined in (6.4) is a bijection from the fundamental paral- 
lelogram in C with vertices 0, AK, 2iK' and AK + 2iK' to the curve T, it brings 
the group product (6.3) down to addition on the torus, f(y) ■ f(z) = f(y + z). 

Because the group is abelian, the subset of V defined as A := {p G Y \ p 2 = e} 
is also a subgroup. Explicitly 

A = {(0, 1), (0, -1), (1/e, -l/e 2 ) £ ^ , (-1/e, l/e 2 )^ } , (6.5a) 
r : (A) = {0,2K,iK',2K + iK'}. (6.5b) 

In what follows a technical, but important role will be played by this subgroup. 
One reason for this is that if p G A, q e A or p • q _1 E A, then the polynomial 
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(6.1) is reducible (in the projective sense, i.e., the leading order term of the 
polynomial is reducible when the unbounded points are considered). In partic- 
ular this means that the Backlund transformation of the equation defined by 
(6.1), with Backlund parameter chosen from A (6.5), reduces to one of the point 
symmetries 

{u — > u, u — > — u, u — > 1/u, u — > — l/u} . (6.6) 

Although its role here is not essential, the three-leg form of Q4, which was 
established in [6] and given for the Jacobi variables in [23] , will be useful. Intro- 
ducing the uniformizing variables a, 8 and £ by writing p = f(a), q = f(/3) and 
u = Vfcsn(^), we have the following identity 

- , _ „ ~. fc 3 \/fcsn(£; fc)sn(a; fc)sn(/3: fc)sn(a — 8: fc) 
Q p , q (u,u,u,«) = — -x 

([sn(£ + a;k)- sn(|; fc)] [an(£ - 8; fc) - sn(£ fc)] [sn(£ - a + /3; A) 
- [sn(£ - a; k) - sn(|; fc)] [sn(£ + /3; fc) - sn(|; fc)] [sn(£ + a - 8; k) 

in which we have introduced 7 given by 

1 - fc 2 sn(a; fc)sn(/3; fc)sn(£; fc)sn(£ + a — 8; k) 
^ 1 — fc 2 sn(a; fc)sn(/3; fc)sn(£; fc)sn(£ — a + 8; fc) ' 
_ 9(g + a; fc)0(g - 8; fc)6(g - a + fc) 
~ 9(£ - a; fc)9(£ + /?; fc)6(£ + a - /?; fc) ' 

The second way of writing 7 in (6.8) which involves the Jacobi 9 function (see 
appendix) leads to the three-leg form, the equality between the two expressions 
for 7 is an elliptic function identity ((A. 5) in the appendix). Assuming the first 
expression for 7 in (6.8) the identity (6.7) may be verified using only the addition 
formula encoded in (6.3). 



- sn(£; fc)] 7 

-sn(f;fc)]), 
(6.7) 

(6.8) 



6.2. Construction of the solution v. The solution for QA given in [4] was found as 
a 'fixed-point' [18,19] or '1-cycle' [17] of the Backlund transformation. That is, a 

solution which is related to itself, v ~ v for some freely chosen Backlund parame- 
ter t = (t, T) e T\ A. In Definition 2 we introduced the extended Backlund trans- 
formation. However for the moment we adhere to [4] and use the un-extended 
version, so that the defining equations for v are 

Qp,t(v,v,v,v) = 0, Q Pit (v,v,v,v) = 0. (6.9) 

These are compatible symmetric biquadratic correspondences, their solution is in 
terms of shifts on an elliptic curve which, whilst being again of Jacobi type, has 
an elliptic modulus different than the modulus fc associated with the equation. 
Thus, in order to write the solution of (6.9) we need to introduce a second elliptic 
modulus, fc* , associated with the solution, v. The new modulus is defined in terms 
of the free parameter t by the relation 

1 1 — T 

k * + TT = 2 ^- ( 6 - 10 ) 
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In addition to the new modulus we also define a new mapping 

f* : z (^y/k^sn(z; cn(z; fc»)dn(z; k*)j , (6-11) 

which is a bijection from the fundamental parallelogram in C with vertices 0, 
4iC, 2iKl, and &K* + 2iK' Sfl to the new elliptic curve r* = r(fc*). Note that 
in (6.3) we gave the rational representation of the group product on Y indepen- 
dently of the modulus k, so the product on T* has the same rational represen- 
tation, and of course f*(y) ■ f*{z) = f*(y + z). 

Now, in order to write the solution of (6.9) explicitly, we also need to introduce 
a relation between T and r*. Specifically 

S:={((p,P),(p*,P*))erxT*\ P l=pf^,P* = ]( P -^^)}. 

(6.12) 

The relation 6 does not establish a bijection between the two curves, in fact a 
generic point on one curve is related to two points on the other. 3 In order to 
write down the solution of (6.9) we introduce two points in r*, 

p* = (p*,P*) = f*(a*), q* = (q*,Q*) = f*{0*) (6.13) 
which are defined in terms of the lattice parameters p and q through S, 

(p,p*).(q.q.) e*. (6.14) 

The solution v is then given by 

v = \fkl sn(£*; fc«), = + na * + mfi*-, (6.15) 

where £q is an arbitrary constant. 

That the function v in (6.15) satisfies (6.9) was demonstrated in [4]. Here our 
main concern is that this function satisfies Qp t q(v, v, v, v) — 0, i.e., is a solution 
of the equation defined by (6.1). This can actually be verified quite easily by 
first checking the following identity: 

Qp, q (u,u,u,u) = ^(—-—)Q -i(u,u,u,u) + ^(— + —)Qp,, q .(u,u,u,u), 

(6.16) 

which holds on any function u = u(n,m), relying only on (6.14). It turns out 
that each term in (6.16) vanishes on the function v given in (6.15). This can be 
seen by using (6.7), where the substitution (u,p,q) — > (v, p*, q J 1 ) leads to the 
identifications (k,a,(3, £) — ► (fc„, a*, ±/?*,£*), whilst from (6.15) £* = and 
so on. (The function v is the 'non-germinating seed' [4] or 'singular' [7] solution 
of the equations Q ±i (u, u, u, u) = 0.) 

To finish this subsection we give the covariant extension of the solution v, 
which is necessary in order to apply the AT-soliton formula (4.4). This is very 
natural because v was constructed as a 1-cycle of the Backlund transformation, 



3 Consideration of the definition (6.12) reveals that (p, p») 6 <5 -4=> (t ■ p -1 ,p«) g 5 <t4> 
(p, pjT 1 ) £ S, therefore 8 establishes a bijection between the sets {{p, p -1 • t} | p 6 T} and 
{{p», p^ 1 } | p* 6 r»}. The group structure on T» turns the latter set here into an abelian 
2-group [22], so <5 endows the former set with this structure as well. 
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in fact a covariantly extended solution emerges if we construct the 1-cycle of the 
extended Backhand transformation (3.2). This amounts to complementing the 
system (6.9) with the equations 

ei i ,t(«,T j t;,«,T i i;) = ) ie{l...N}, (6.17) 

where we have introduced the Backhand parameters [j = (k,Li) e T, i e 
{1 . . . N} and the shifts Ti . . .Jn- The function v in (6.15) may be extended to 
solve the system (6.17) explicitly by introducing the points I* = (l*,L*) = f*(A*), 
ie {I...N} such that 

(k,l*)eS, ie{l...N}, (6.18) 
and complementing (6.15) with the (evidently compatible) equations 

T^=^+A|, ie{l...N}. (6.19) 



6.3. Integration of the equations for <px . . .(j>N- In the previous subsection we 
have recalled from [4] an explicit solution, v in (6.15), for the equation Q4, and 
given its natural covariant extension through (6.19). To apply the main theorem 
in Section 4 it remains to give the functions <pi . . . (f>N which we construct from 
v by integration of the equations (4.2). It turns out that with Q as in (6.1) and 
v as in (6.15), the system (4.2) for fa may be reduced by the substitution 

9(^ + 2A*;M 

* = Pl ( } 

to the following autonomous system for the new variable pi A : 

p,i,+ P (*je(A 1 '+a,;t)''" Pl \qJi + ql*J + f3*;K) PU 



(6.21a) 



We can integrate the equations (6.21a), so that 

ff p*k -pl* \ Q{\* - a,;k,) \ n ff qjj - gl* \ e{\*-p*;K) 
Pl Pl '°{{pJi+pl*)e(\*+a*;k*)) \\qj l + ql*)e{\*+p*;k*) / 

(6.22) 

for i e {1 . . . N}, where pi, ■ ■ ■ Pn,o ar e independent of n and m. Thus for each 
i£ {1 . . . N} the function </>j is given by (6.20) in terms of (6.15), (6.19) and 
the new function pi (6.22), (6.21b). 

It remains to satisfy the reader that the substitution (6.20) reduces the system 
(4.2) to the autonomous system (6.21). Let us focus on the first equation of 
(4.2a), the other equations present in (4.2) are similar. We begin by using the 



4 The p variable introduced here is different from the one defined in [4] . 
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identity (6.16) which leads to the following reformulation of the first equation in 
(4.2a) ' 

7 Qp,h{^iV,^iV,l1v,v) 

<Pi = -7: Tt =r2^T <Pi 

Qp,u{TiV,TiV,v,Tfv) 

pJi-pl*\ Q Pt , n -i(T t v 1 T l v 1 T 2 l v 1 v) ' '- ,J 



p*k+pl*J Q p ,,i*(TiV,TiV,v,T?v) 

That <2 P „,[* (TiV, TiV, Tfv, v) = and Qn^ir-i {JiV, TiV, v, T?{5) = has also 
been exploited here, these are verifiable using the expansion (6.7). We then 
propose the following string of equalities 

Q Pt ,[«0>,T*'M,T!w) 

1 - fcgsn(A*;fc„)sn(A* + a*; fc»)sn(£„ + A*; fc„)sn(£„ + a* + A*;fc») 
1 - £;2sn(A*;fc*)sn(A* - a »; fc*)sn(£* + A*; fc st )sn(^» + a* + A*;fc*)' 

_ 9(A* - a,; fc*)9(6. + 2A* + a,; fc,)6(^; fc,) 

~ 6(A* + a*; fc»)©(£* + a*; fc*)6(^ + 2A*; K) ' 

_ 6(A* -a.ifcQefe +2A?;fc„)e(g.;fc.) 

~ 6(A* +a»;fc»)e(6;fc»)e(^ + 2A*;fc*)' 

(6.24) 

The first equality here can be verified using (6.7) (modulo some manipulation). 
The second follows from the first by an addition formula ((A. 5) listed in the 
appendix) . The third equality is immediate from the second because = + a* 
by (6.15). Finally, combining (6.24) with (6.23) motivates the substitution (6.20) 
and clearly results in the first equation of (6.21a). 



6.4- The N -soliton solution. In Sections 6.2 and 6.3 we have given the ingredi- 
ents for the ./V-soliton formula (4.4), that is a solution v, its covariant extension, 
and the functions <pi . . . (j)^. To write down the iV-soliton solution constructed 
from these ingredients we consider first the function / defined in (4.14) which 
appears in the denominator of (4.4). Substituting (6.20) and (6.15) into (4.14) 
whilst bearing in mind the expansion (4.15), and subsequently using (6.21b) and 
(6.19) to write the result in terms of the un-shifted functions and pi . . . pn, 
we find 



y^ e(6»+2A*;fc*) p , y> +2\* + 2\*;k*) ^ x 

i=i ©(£*;&*) »=ij=»+i @(£*;^*) 



^ A A e(£* + 2A*+2A* + 2A*;A;*) 
"EE E~ K7>-n —PiPiP^lXf K X^ + ..., (6.25) 

where we have introduced the constants 

f lJU - Ijlt \ 9(A* - A*; fc.) _ ( 1*1* - hi, \ H(A* - Xjjk*) 

11 ■ [ri, ■ / ; /;A-)iA; • a; : /,m ■ /,/JiisA; ( ' 



A*; fc*) 
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for each i, j e {1 . . .N}. These constants appear because TjP% = XijPi as can 
be seen by comparing (6.26) and (6.21b). The equality between the alternative 
forms we have given for Xij in (6.26) can be verified using (A. 2) together with 
the addition formula encoded in (6.3) and the relations (6.18). A similar consid- 
eration of the numerator in the A-soliton formula (4.4), whilst bearing in mind 
(A. 2), makes it natural to introduce a new function g, 

-tt t m ' + ^;^ k ' } ^^ + .... „ 27 , 

i=l j=i+l K=j+1 *' 

The A-soliton solution for the equation defined by (6.1) may then be expressed 
as 

u {N) = vj; (6.28) 
with v, f and g as in (6.15), (6.25) and (6.27), explicitly for N e {0, 1, 2} this is 

(o) _ H(C*i fc*) (i) _ H(^; fc^) — H(^ + 2A 1 ; k*)p\ p) _ 
_ e(£*;fc*)' U ~ e(6;fe.)-6(^+2AI;A;,)pi' " 
n(^]K)-R(^ + 2Xl;K)p 1 -B{^ + 2X*-,K)p 2 +li(^+2X* 1 + 2X* 2 -,K)p 1 p 2 X^ 
9(6*; fc„)-9(e* + 2A|; K) Pl -<d(^ + 2X* 2 ; K)p 2 + Q(£_ r + 2X* 1 + 2X* 2: K) Pl p 2 Xf 2 ' 

The functions £* and p\ . . . pn appearing in the solution are given in (6.15) and 
(6.22). The modulus k and the constants p = (p, P) and q = (q,Q) taken from 
r = r(fc) defined in (6.2) are fixed (they appear in the equation), whereas the 
constants t = (t,T) e V \ A and U = (U,Li) = f(A;) € T, p ija £ CU {oo} for 
i e {1 . . .N} may all be chosen freely. The modulus fc* is determined from t 
by (6.10), whilst the parameters p* = (p*,P*) = f*(a»), q* = {q*,Q*) = f*(/3*) 
and [* = (/*,£*) = f*(A*) fori e {1...N} which lie in r* = T(k*) are all 
determined from their counterparts in T by the relation 5 defined in (6.12), 

(p,p»),(q,q»),([i,[i)-..([Ar,nv) G 

To ensure that the A-soliton solution (6.28) has distinct functional depen- 
dence on each of the N constants of integration pi o . . . pn,o appearing in (6.22), 
and so constitutes a true A-soliton solution, we should take some care in how 
we choose the parameters li . . . In. First of all we should choose ti . . . In € 
r \ (A U A ■ t), this restriction ensures that, for each ie {1 . . . A}, the Backhand 
equations associated with parameter [, are not reducible, and also that the covari- 
ant extension provides two distinct solutions of the Backlund equations applied 
to the solution v, i.e., TiV ^ T 1 ?;. Furthermore, because we construct higher 
soliton solutions by superposition, we should also stipulate that U ■ IJ 1 A for all 
distinct i,j e {1 . . .N}. Let us remark though, that soliton solutions for which 
[, £ A ■ t for some i € {1...N} or U ■ lj l e A for some distinct i, j G {1 . . . N}, 
are constructable by quadrature 5 , or recoverable from the solution given here by 
taking a subtle limit. 



5 Although the main results here rely on there being two, a single particular solution is in 
fact sufficient to subsequently solve a Riccati equation by quadrature. 
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6.5. Special lattice directions. Rather than avoiding them altogether, we will 
now make particular use of the elements in A. It is quite straightforward to see 
that rnr» = A, so A is a subgroup of both T and T* 6 . It turns out that the 
relation 5 (6.12) is compatible with this shared subgroup structure on the two 
curves, p e A => (p,p) € S. We will focus on the element e e A as well as one 
other element which for convenience we denote by [ i , 

Ii := (1/e, -l/e 2 )_ = f(itf') = WK)- (6-29) 

Note that for a generic point p = (p, P) £ T, p ■ [i = (1/p, —P/p 2 ). We now 
introduce two elements of S, 7 

(t,e), (t-[i,Ii) G S (6.30) 

and to each we associate a new special lattice direction into which we extend 
and pi ... pat, denoting the shifts by To and Ti. Specifically (using the properties 
listed in the appendix) 

T £*=£*, T ^* = ^* + iK' t , 
Tqv = v, Tiv = 1/v, 

To4>j =~<f>j, Tifa = - (sn(^* +2A^;fc*)/sn(^*;fc*)) <j>j, 

for j e {1 . . . A^} (again, note that in (6.31) i denotes the imaginary unit). Our 
primary interest in these special lattice directions is that they facilitate the 
construction of a new operator 

X := wT Ti (6.32) 

which has the two properties 

XI = v, [Bj,X]=0, je{l...N}. (6.33) 

Using these properties in (4.4) we are able to establish a rather simple relation- 
ship between the functions / and g appearing in the A^-soliton solution (6.28), 
namely 

<? = T Ti/, (6.34) 

which can also be verified directly using (6.31) and the expressions for / and g 
(6.25) and (6.27). 



6 This is closely related to a characterisation of the deformation r — ► T* which was estab- 
lished in [24]. 

7 Actually we could choose the points (e, e), (li , G S here, which at first seems more 

2 2 

natural, however the subsequent calculations require more effort. 



20 



J Atkinson and F W Nijhoff 



7. Discussion 

We have described a technique which enables solution of the Backlund iteration 
scheme for a class of integrable lattice equations. To express the solution after 
TV applications of the Backlund transformation we require that the covariant 
extension of the initial solution be known. This is a natural concept emerging 
from the feature of the lattice equations that lattice and Backlund directions are 
distinguished only by the difference of a parameter, and therefore the technique 
can be said to have a manifestly discrete origin. The relative mildness of the 
covariant-extendibility requirement means that this approach usefully separates 
two hitherto intertwined aspects of the soliton-type solutions; the choice of seed 
solution on the one hand, and the solution of the Backlund iteration on the 
other. 

The iV-soliton formula we have discovered holds for all equations listed by 
Adlcr, Bobenko and Suris (ABS) in [6] and is independent of the particular form 
of the polynomial defining the equation. It can be said to have Hirota form, 
but the technique also introduces some new features, in particular there is a 
commuting family of second-order linear difference operators naturally associ- 
ated with the solution, and these operators essentially factorise the Hirota-type 
polynomial. 

An explicit Af-soliton solution of the primary model listed in [6], Q4, has been 
given as a particular instance of the A^-soliton formula. We expect this solution 
can be written (and directly verified) in a Cauchy-matrix form generalising re- 
sults in [5], in fact an elliptic Cauchy-matrix approach has been developed by 
the authors alongside the present work [25]. Direct verification of this solution 
through bilincarization and a Casorati determinant form as in [26] , which would 
generalise results of [9], is also likely to be possible. 
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A. Jacobi elliptic and theta functions 

In Section 6 the equation considered and its solution are given in terms of the 
Jacobi elliptic and theta functions. Definitions of these functions we adhere to 
can be found in Chapter 5 of [21]. For the convenience of the reader this appendix 
recalls some of the properties of these functions. 

By convention the elliptic modulus is denoted by k and the half-periods of 
sn(z) = sn(z; k) by 2K and iK' . Increments in the argument by the half-periods 
result in the following behaviour 



The function sn(z) may be written in terms of the Jacobi theta functions H(z) = 
R(z;k) and Q(z) = Q(z;k), 



sn(z + 2K) = -sn(z), sn(z + iK ') = 



1 



(A.l) 



ksn(z) 




H(z) 



(A.2) 



e( z y 



(A.5) 
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In turn these functions satisfy the basic relations 

R(z + 2K) = H(-z) = -H(z), 
Q(z + 2K) = Q(-z) = Q(z), 

H(z + iK') = ie-^ 2z+lK ')/ 4K e(z), (A ' 3) 

@(z + iK') = ie-*< 2z+tK 'V iK K(z), 

where i is the imaginary unit. We also note the important addition formula 

U(x + y)R(x - y)K(z + w)K(z - w) 
+H(x + z)~R(x - z)~R(w + y)R(w - y) (A.4) 
+H(x + w)R(x - w)U(y + z)U(y - z) = 0, 

from which others can be derived. One identity in particular is used in the text: 

1 - fc 2 sn(a)sn(6)sn(c)sn(c + a - b) _ 0(c + a)6(c - b)Q(c -a + b) 
1 - /c 2 sn(a)sn(6)sn(c)sn(c -a + b) ~ 9(c - a)6(c + b)Q(c + a - b) ' 

To derive (A.5) we choose 

w = c-(a-b)/2, x = c+(a-b)/2, 

y=(b-a)/2, z = (a + b)/2 + iK' 

and 

w = c-(a- b)/2 + iK', x = c+(a- b)/2 + iK', 

y=(b-a)/2, z = {a + b)/2 

in (A.4), leading respectively to 

9(c + a)6(c - 6)H(c - a + 6) - 6(c - a)6(c + 6)H(c + a - b) 

= H(2c)H(6-a)e(a)e(6)/H(c) (A.6) 

and 

6(c + a)6(c - 6)6(c - a + 6) - 6(c - a)6(c + 6)6(c + a - 6) 

= H(2c)H(6-a)H(a)H(6)/9(c), (A.7) 

which we have simplified using (A. 3). Elimination of H(2c)H(6 — a) between (A.6) 
and (A.7) followed by the use of H(z) = sn(z)9(z) from (A. 2) to remove the 
remaining occurrences of H results in (A.5). 
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